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ABSTRACT 

At high luminosities black hole binaries show spectra with a strong disc component 
accompanied by an equally strong tail where at least some of the electrons are non- 
thermal. We reanalyze the simultaneous ASCA-RXTE-OSSE data from the 1998 out¬ 
burst of XTE J1550-564, which span 0.7-1000 keV and remain the best data available 
of a black hole binary in this state. We reassess the importance of electron-positron 
pair production using a realistically high value of the source compactness for the first 
time. The lack of an observable annihilation line together with the observed y-ray 
flux beyond 511 keV constrains the maximum electron Lorentz factor to be <10 and 
the slope of the injected electrons to <2.5. We also use the fast (10-50 Hz) variability 
spectrum to constrain the spatial dependence of the electron heating and acceleration. 
We find that the spectrum of the fast variability is consistent with being fully thermal, 
so that the observed non-thermal emission is produced from further out in the flow. 

Key words: Accretion, accretion discs - X-rays: binaries - X-rays: individual (XTE 
J1550-564) 


1 INTRODUCTION 

The so-called very high state is perhaps the most complex 
and least understood spectral state of black hole binaries. It 
is a highly luminous state characterized by a strong black- 
body component together with a very strong rather steep 
tail, typically starting at the peak of the disc component 
which can make these two components hard to separate from 
each other. The question of whether the disc in the very high 
state is truncated or not is a matter of debate. It is also 
not clear whether the strong signs of non-thermal Comp¬ 
ton emission is the result of scattering off one single hybrid 
electron distribution, or if the thermal and non-thermal elec¬ 
trons come from different spatially separated regions of the 
flow. A popular idea of the possible configuration is that the 
thermal electrons are located close to the black hole in a 
region which may be the remains of the hot inner flow seen 
in the low/hard state, while the non-thermal electrons are 
located further out in the flow, perhaps accelerated by mag¬ 
netic flares above the disc. To distinguish between a single 
hybrid distribution and a thermal plus a non-thermal (or to 
be correct, one with purely non-thermal injection, see dis- 
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cussion in Section 3) is not trivial but is of important for 
understanding the heating/acceleration processes involved. 

Further, the persistent non detection of any annihilation 
lines from Galactic black holes questions the importance of 
pair production in these sources, and is often used for as¬ 
suming a strong upper limit on their compactness. While 
this may be justifiable in less luminous states, it is difficult 
to motivate in the very high state. 

1.1 XTE J1550-564 

XTE J1550-564 is a transient low-mass X-ray binary. It was 
discovered during a n outburst i n 1998 by the RXTE All- 
Sky Monitor (ASM; ISmith|[T998l f and t he CGRO Burst an d 
Transient Source Experiment fBATSE: I Wilson et al.|[l9^ b 
The companion is a main sequence star of type G-K and the 
mass of the compact object has been estimated to 8.4-11.2 
M© (Orosz02), and it is thus considered to be a black hole. 
The source displays all the canonical spectral states typically 
associated with black hole binaries. The X-ray spectrum and 
its evolution has been studied extensively since its discovery. 

In this paper we reanalyse the simultaneous ASGA- 
RXTE-OSSE data from the 1998 outburst which spans 0.7- 
1000 keV and remain the best broadband data of an X-ray 
binary in the very high state. Previous modeling of the spec- 
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trum has not been able to utilise it to the fullest extent due 
to issues with the ASCA-RXTE cross-calibration and lack 
of good ionized reflection models. For the hrst time we now 
ht the entire broadband 0.7-1000 keV spectrum simultane¬ 
ously. 

We use additional information from the fast (10-50 Hz) 
variability in the form of frequency resolved spectroscopy. 
This is a very powerful tool to help breaking the inher¬ 
ent degeneracy of spectral modelli n g, use d e.g. in the pi¬ 
oneering work by IChurazov et al.l (|200lh . The technique 
has h owever not been used much for t h is pu rpose. Re¬ 
cently, lAxelsson. Hjalmarsdotter Donel (|2013h , hereafter 
AHD13, studied the PCA spectrum of XTE J1550-564 dur¬ 
ing the rise of the strong flare in 1998 and compared the 
spectrum of the continuum to that of the fast (10-50 Hz) 
variability. The difference that was found between the time 
averaged spectrum and the spectrum of the fast variability 
strongly supports spatially inhomogeneous comptonization 
models where (at least) two electron distributions seem to 
be present at different distances from the black hole. In this 
paper we develop this idea further and also attempt to sep¬ 
arate the thermal and non-thermal distributions. 

The broadband data give us a good estimate of the 
source luminosity and compactness, which is high. This 
should lead to a non-negligible production of electron- 
positron pairs. The lack of an annihilation line in the OSSE 
data allows us to derive limits on pair production in this 
source that may be applicable to luminous black hole bina¬ 
ries in general. 

2 DATA 

2.1 PCA data continuum and frequency resolved 
spectra 

We use the data from the plateau phase, approximately two 
weeks after the initial flare. For the PCA continuum data we 
extracted Standard2 spectra for all 8 observations between 
September 23 and Oc tober 6, applying s t andar d selection 
criteria. As shown in iGierliiiski &: Don3 (|2003h . hereafter 
GD03 (their hg. 2), the position of the source in a PCA 
colour-colour diagram during this time corresponds to the 
very high spectral state. ( This state was also la belled the 
extreme very high state bv iKubota Don3l2004l . hereafter 
KD04, to separate it from a more disc dominated very high 
state also displayed by this source.) These were then co¬ 
added to create a total PCA spectrum with a total exposure 
time of 31 ks. A systematic error of 1 per cent was added 
to each bin. The energy band used for spectral modelling is 
3-20 keV. To extract the spectrum of the fast variab il ity we 
follo w the procedure d escrib ed in iRevnivtsev et aP (|l999l l 
and IRevnivtsev et al ] (l200lh . We extract a light curve for 
each available channel for each observation and construct a 
power density spectrum (PDS). The relative contribution 
of each channel was determined by integrating the PDS 
above 10 Hz, and used to construct an energy spectrum of 
the rapid variability. As with the continuum spectra, re¬ 
sults from all observations were then co-added into a single 
spectrum, and a systematic error of one per cent added to 
each bin. We choose the frequency range 10-50 Hz in or¬ 
der to avoid contamination from the observed QPO at < 4 


Hz and its harmonic (see further discus sion in AHD13 and 
lAxelsson. Done Hialmarsdotterll2014l . 


2.2 HEXTE and OSSE data 

We use the HEXTE data, averaged over the same period, 
from cluster 0 only. DG03 found the HEXTE data to be 
consistent with remaining constant throughout the period. 

The OSSE data are from viewing period 729.5, from 
September 25 to October 6, 1998. We use the same high- 
level product spectrum as in DG03, extracted for detectors 

1.2.3 and 4 for the whole viewing period. The OSSE and 
HEXTE data are in very good agreement in the overlapping 
region between 50 and ~120 keV above which the HEXTE 
data start to slightly flatten out. The agreement is however 
still within the HEXTE errors up to 200 keV. The OSSE 
data show no evidence of a high-energy cutoff up to at least 
1000 keV. There is no sign of any annihilation line at 511 
keV in the data. 

Unfortunately the low sensitivity of HEXTE and low 
time resolution with retained spectral resolution of the 
OSSE data do not allow us to extend the fourier-spectrum 
to higher energies. 

2.3 ASCA data 

The ASCA data used here are from a 24 ksec pointing on 
September 24, 1998 and is the same data as used in GD03, 
extending from 0.7-10 keV. In GD03 the ASCA data was 
treated separately to avoid previous cross-calibration issues 
with RXTE. These cross-calibration issues were hxed in 
HEASOFT 5.2 (and in KD04 this ASCA data were mod¬ 
elled together with some of the PCA and HEXTE pointings 
includ ed in our co-added sp ectrum). The data were also used 
e.g. bv ISteiner et al.l (|201lh (see discussion in Section 5). 

3 SPECTRAL MODEL 
3.1 Eqpair 

We use the comptonization code EQPAIR ICoppil (|2000l b The 
key-concept in this model is the compactness, defined as: 


where L is the luminosity or power, R the size of the comp- 
tonizing region, ctt the Thomson cross section and rrie the 
electron mass. 

The hard compactness, /h, is determined by Lh, the lu¬ 
minosity or power supplied to the electrons in the comp- 
tonizing region. This power can be supplied either as ther¬ 
mal heating Lth, or as energetic electrons (and/or positrons) 
being injected into the source with a fixed non-thermal en¬ 
ergy spectrum Q(7) that can be either mono-energetic or 
a power law with index Vinj extending from Lorentz factor 
7mzn ~ 1-2 to 7max ~ 3-1000. The total hard compactness 
is thus /h=hh+Uth- The ratio Rth/R gives the relative im¬ 
portance of non-thermal scattering and if this parameter is 
> 0 we have a hybrid electron distribution with a low-energy 
Maxwellian and a high-energy power law part. 
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The soft compactness, Is, is in turn determined by Lg, 
the luminosity of soft photons available for Compton up- 
scattering and thus cooling of the hot electrons. If electron 
heating/acceleration dominates over cooling we get a hard 
spectrum and if cooling by soft photons dominates the re¬ 
sulting spectrum is soft, and the main parameter of the 
model determining the overall spectral state is thus the ratio 
Ih/ls- The absolute value of the soft compactness Is is also a 
fit parameter. This parameter not only measures the effec¬ 
tiveness of Compton cooling, compared to Coulomb cooling, 
but also the optical depth for pair-production in the source. 
Note that in EQPAIR Is is not coupled to the normalization. 
In practice, the fit is therefore insensitive to this parameter 
(except regarding the strength of the annihilation line, see 
further discussion in the next section) and it is therefore of¬ 
ten kept frozen to some fiducial value in the fitting process. 
The exact spectral shape for a given amount of seed photons 
is determined by a combination of the electron temperature 
and optical depth of the hot flow. The total optical depth 
tt may in addition to the electron optical depth Tp, which 
is a fit parameter, also include a part from pairs. Together 
with the electron temperature kTe, tt is calculated by the 
model self consistently. The model does not automatically 
assume pair balance but pairs that are injected or created 
in the plasma annihilate away once they have cooled. 

The model assumes a spherical geometry with the pho¬ 
tons injected homogeneously throughout a spherical comp- 
tonizing cloud. As input source of soft photons EQPAIR in¬ 
cludes a soft photon source modelled as a disc or single¬ 
temperature blackbody. Since the seed photon temperature 
is mainly decided by the maximum temperature disc pho¬ 
tons, we use the blackbody option for simplicity. The part 
of the photons reaching us un-scattered is normally given in 
output by the model without the addition of any extra disc 
component. We have here however modified EQPAIR for use 
with a convolution reflection model and removed the black¬ 
body from the output spectrum to avoid having the disc 
reflected in itself. We therefore model the direct disc emis¬ 
sion with a separate disc component using DISKBB. Following 
KD04, who investigated the temperature-luminosity correla¬ 
tion in spectra based on the same ASCA-PCA-HEXTE data 
as we use here, we assume that the disc is slightly truncated 
and using other more sophisticated disc mo dels including 
relati vistic effects thus makes no sense (but see lSteiner et al.1 
120111 for a different view). The maximum temperature of the 
disc is allowed to differ from that of the input photons in 
EQPAIR since in reality we may see a somewhat cooler part 
of the disc in direct emission than that providing the seed 
photons. 

3.2 RFXCONV and the iron-line region 

Here, as in AHD13, we use the relatively new model RFX¬ 
CONV to model reflection of the hard X-rays off the ac¬ 
cretion disc. The RFXCONV model is a convolution model 
that combines t he ta ble models for the reflected spectra of 
IRoss EabianI (|2005l j from a constant density ionized disc, 
with the IREFLECT con volution variant of the PEXRIV Comp¬ 
ton reflection code bv iMagdziarz ZdziarskJ (|l995l b It in¬ 
cludes a self-consistently calculated iron line, convolved with 
any continuum shape. The main parameters are the relative 
amplitude of the reflected component, R, the inclination and 


the ionization parameter of the reprocessing matte r, Eor 
the inclination we use i — 70° (jOrosz et al.ll2002l b We al¬ 
low for relativistic smearing of the reflection features using 
KDBLUR. We include photoelectric absorption in the form of 
TBABS. In AHD13 the fits to the PCA data were found to 
be significantly improved by the addition of a narrow ab¬ 
sorption line (gabs) at ~ 6.8 keV. Observations of similar 
strong absorption lines have been reported in several high 
inclination Galactic binary systems and are believed to rep¬ 
resent iron K resonance at 6.7 and 7.0 keV from outflowing 
material fe.g. JUeda et al.iri998l : iLee et ahllioo^ : lUeda et al.l 
l2004h . We thus include a line as GABS with central energy 
between 6.4 and 7.1 keV. 

Our initial single-component hybrid 

comptonization model of the form GON- 
stant*tbabs*gabs(diskbb+kdblur*rfxgonv*eqpair. 
Compared to the similar HYB model in GD03 we use a 
broader energy range including the ASCA data, a better 
reflection model RFXGONV and different assumptions about 
the compactness, see below. 


4 RESULTS 

4.1 Single-component hybrid comptonization 
model for the 0.7—1000 keV continuum 

Eigure la shows the best fit total model of the form GON- 
STANT*TBABS*GABS(DISGBB+KDBLUR*RFXGONV*EQPAIR)to 
the 0.7-1000 keV continuum data of XTE J1550-564. The 
parameters are listed in Table 1, left column (model ONE- 
GOMPONENT HYBRID) and the components of the model are 
plotted in Eig. lb. 

The spectrum contains a strong disc component. The 
best-fit inner disc temperature of the unscattered disc com¬ 
ponent is 0.56 keV, in agreement with or slightly higher than 
in DG03, KD04 and AHD13. We use the disc temperature 
as input for the seed photon temperature but initially allow 
them to be different. The best value for the seed tempera¬ 
ture however stays the same as that of the inner disc and we 
thus set them equal in order to reduce the number of free 
parameters. 

We also see a very strong Comptonized component, typ¬ 
ical for the very high state. Our best fit value of /h/4~ 1 
indicates equal power in heating of the electrons as in cool¬ 
ing them off soft photons. The electron distribution is hybrid 
with ~ 50% of the power to the electrons supplied as accel¬ 
eration as opposed to thermal heating. The temperature of 
the low-end Maxwellian electrons in our best fit is kT^ — 7.0 
keV. The optical depth is tt = 3.9, and includes a small 
contribution (r = 0.1) from pairs. We stress that in a soft 
spectrum with hybrid electrons present the rollover from the 
low energy Maxwellian distribution falls in the same energy 
range, 3A;Te ~20-30 keV, as the peak of the reflection bump 
and is therefore difficult to determine exactly. The exact 
values of the electron temperature and optical depth may 
therefore be somewhat degenerate within given errors. 

In our best fit, reflection is moderately high at 0.65 and 
comparable to values found in AHD13 and GD03. Reflection 
is highly ionized with best fit value log ^=2.7. The inclusion 
of a broad absorption line with a = 0.40 centered at 6.8 keV 
significantly improves the fit. We are not sure as to whether 
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Figure 1. One-component hybrid Comptonization model to the 0.7-1000 keV data of XTEJ1550-564. Left panel: Data and model 
including residuals. ASCA data 0.7-10 keV in black, PC A 3-20 keV in red, HEXTE 20-200 keV in green and OSSE 50-1000 keV in 
blue. Right panel: Components of the model. The unscattered blackbody in magenta long dashes, Comptonization from a hybrid electron 
distribution in blue short dashes and Compton reflection as red dots. 


this line is real or if it is an artefact compensating for the 
strength of the fluorescent line as compared to continuum 
reflection in the RFXCONV model. We stress that in RFXCONV 
where the iron line is calculated self-consistently, it can not 
be arbitrarily moved in energy or strength (as when using 
e.g. IREFLECT + gauss). 

We include KDBLUR but find that the fit is insensitive to 
the value of the inner radius together with ionized reflection. 
We find that values lower than = 170 are ruled out by 
the data and thus freeze the inner disc radius to this value 
rather than to a best-fit value even if tends to decrease 
for even larger values. We do not believe this value to repre¬ 
sent a true lower limit on the inner disc radius. But it shows 
that a small inner radius with signs of relativistic effects on 
the reflection features is not compatible with the data if we 
assume ionized reflection. We find that the only way our 
data is compatible with small inner disc radius is if we as¬ 
sume zero ionization of the reflector. This does not only give 
a worse fit statistically, we also have good reason to believe 
that the reflection features should show signs of ionization 
since the disc is hot as well as strongly irradiated. The ion¬ 
ization parameter is rather a representation of an average 
than the true ionization state of the reflector as the reflection 
spectrum is a composite of many different parts of the verti¬ 
cal structure of the disc. Some contribution comes from the 
almost completely ionized surface while some comes from 
the almost completely neutral material further inside the 
disc, but very little reflection is expected from intermedi¬ 
ate ionization states ([Navakshin. Kazanas Kallman|[200g : 
IPone Sz Navakshinll2QQ7l h 

4 . 1.1 Compactness and annihilation lines 

There is no sign of an annihilation line at 511 keV in the 
data of XTE J1550-564. The lack of observed annihilation 
features in GBH data is often used as motivation for con¬ 
straining Is < 10. While this may be .justifiable in l ess lu - 
minous sources, like e.g. Cyg X-1 (|Gierlinski et al.lll997li . 
it is certainly not so in the more luminous sources and/or 
states. Our model gives a soft luminosity of ~ 1.7 x 10^^ 


erg s“^. In the very high state, we expect a fairly large 
overlapping region of the hot flow and the disc. But a soft 
compactness of ~ 10 would imply a comptonizing region ex¬ 
tending beyond 500 Rg which is not easily reconciled with 
our current picture of the accretion geometry. More reason¬ 
able radii of the extension of the hot flow, 10-30 Rg, give 
a range of /s~200-700. Since the modelling is insensitive 
to this parameter (in all but the size of the annihilation 
line) we thus keep Is frozen to the more realistic value of 
500. Such a high compactness should lead to copious pair 
production given that enough high energy gamma rays are 
present. Using a realistic value for the compactness based 
on the luminosity together with good-quality high-energy 
data, we are thus in a position to constrain the injection 
spectrum of high energy electrons, determined by a combi¬ 
nation of the power law index of injected electrons Finj and 
their maximum gamma factor Xmax- The parameter value 
for the maximum energy of the injected electrons Xmax is 
insensitive to the fitting process and should be kept frozen. 
To determine allowed combinations of T inj and Xmax we fix 
Xmax to a range of different values between 3 and 1000 and 
calculate error ranges for Finj for each Xmax- Fig. 2 shows 
examples of different combinations of T inj and Xmax all giv¬ 
ing acceptable fits (Ax^/^^ < 1) to the broadband data, but 
with a preference for values around Xmax ~ 10 and Finj 
between 2.10-2.50 with best-fit value 2.27. 

In addition, the OSSE data give an upper limit for the 
equivalent width EW of the line of 70 keV (at 3 cr confi¬ 
dence). Fig. 3 shows best-fit models for different combina¬ 
tions of Finj and Xmax from Fig. 2. We find that all models 
with 7ma®>10 produce an annihilation line with EW>70 
keV and can thus be ruled out based on the non-detection 
on basis of the OSSE data alone. We thus keep Xmax frozen 
to 10 (since this parameter is rather insensitive to the fitting 
process as discussed above). 

4.2 The spectrum of the fast variability 

We now fit the spectrum of the 10-50 Hz variability for the 
co-added PGA observations. Since the data start only at 3 
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MODELS 


ONE-COMPONENT 

TWO-COMPONENT 

TWO-COMPONENT 



HYBRID 

THERMAL HYBRID 

THERMAL FULL NON-TH INJECTION 

Ah 

lO^^cm ^ 

0-O0_o.02 

0-O0_o.02 

0 

’0-'J-‘--0.02 

Eabs.line 

keV 

6 77+0-07 

o- ‘ ‘ -0.06 

c 77+0.08 
‘ ‘ -0.05 

- 

*^abs.line 

keV 


'^•^'^-0.12 

- 

Line depthabs.line 

keV 

n 10+0.07 

o-J-O_0.06 

0 14+0-07 

- 

Pdisc 

keV 

^•OO_0.03 

U.bD-0 04 

0 >5fi+0-01 

U.UO_o 02 

Pseed,var 

keV 

- 

0.56/ 

0.56/ 

^h/^s,var 


- 

111 + 0.02 
^•-^^-0.04 

1 1 1 + 0.02 
-^•^-^-0.04 

kT® 

e,var 

keV 

- 

5.46 

5.46 

' p,var 


- 

zLS+O-^O 

6-48tS;J^ 

T,var 

Reffvar 


- 

6.48 

6.48 


- 



Tseed 

keV 

0.56/ 

0.56/ 

0.56/ 

^hAs 


0.981^2 

Q 97+0-04 
-0.04 

n sq+O-08 

Is 


500f 

500f 

500f 

Inth/lh 


0-481^3 

’0-O"±_o.05 

l.Of 

kT« 

keV 

7.00 

6.96 

5.49 

Tp 


' ' -0.17 

3 ^Q+O-^^ 

d-OH-o 47 

9 o/:; + 0.59 

Z.dD-o 24 



3.93 

3.77 

2.51 

F- • 

^ inj 


r\ r\iy-\-0.23 
-0.17 

2 4c;+0-55 

^•^0_o.23 

q 7 [;: + 0.23 

-J- ' 0-0 43 

Reff 


U.DO_q 22 

U.DD_q 24 

1 - 0 - 0.21 

log ^ 


9 7n+0.09 

'O_o.l4 

9 70 + 0.17 
' '^-0.09 

0 /< 7 + 0.06 
-0.04 

^in 

rg 

ITOf** 

170f 

170f 

X^/dof 


245/290 

240/291 

327/295 


'^Calculated self-consistently by the model. 
^Lower limit 
^Parameter fixed. 


Table 1. Fit results of three different models for the continuum as described in Sections 5.1 and 5.3 . Parameters with subscript ’var’ 
represent best-fit values for the thermal model for the spectrum of the fast variability as described in Section 4.2. Errors for these 
parameters are from the fit to the variability spectrum only. They were frozen to their best-fit values when determining the other 
parameters in the two-component models. The unabsorbed bolometric luminosity is Ltot = 3.7 x 10^® for all models. 


10 50 100 500 1000 



Figure 2. Contour plot for Pinj versus 7max from the one- 
component model fit to the data of XTE. Shaded regions show 
contours corresponding to (from darker to lighter colour) Icr, 3a, 
5(7 and 7a. 

keV, we keep the disc temperature frozen to the best fit value 
from the continuum fit based on the ASCA data. There is 
no sign of a blackbody or disc component in the spectrum of 


the fast variability. This agrees with the findings in AHD13 
and confirms the well known result that the fast variability 
does not come from the accretion disc in soft states (e.g. 
IChurazov et aklBoOlh . We thus keep the disc normalization 
frozen to 0. 

Even after removing the disc component the continuum 
model still does not give a good fit to the spectrum of the 
variability. It is clear that the two spectra have a different 
shape in the overlapping energy region (3-20 keV). The spec¬ 
trum of the variability is harder with a higher Ih/ls, meaning 
that the variable part of the flow sees less soft photons than 
the overall accretion flow. The optical depth is also higher. 
This result agrees with AHD13 who showed that the spec¬ 
trum of the fast variability had a different shape, always 
harder than the continuum, in all PC A observations in the 
rise of the strong flare in 1998. 

Since our variability spectrum continues up to 20 keV 
only, we can thus not constrain the contribution of non- 
thermal comptonization or get a very good idea of the elec¬ 
tron temperature from these data. We then have two obvious 
choices regarding the non-thermal electron distribution in 
the spectrum of the fast variability. One would be to freeze 
the ratio of non-thermal to thermal contribution and the 
parameters of the injection spectrum to those of the contin¬ 
uum, thus assuming the same non-thermal electron distri- 
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Figure 4. Spectrum of the 50-100 Hz variability in XTEJ1550-564. Left panel: The fourier resolved PC A data and best-fit (fully thermal) 
model including residuals. Right panel: The model components are shown as dashed lines for the comptonized spectrum and dotted lines 
for its reflection. 



Figure 3. Annihilation lines. Example best-fit models for dif¬ 
ferent combinations of Vinj and ^max- From bottom to top 
7max=10, 20, 50, 100 and 1000. Note that all these models give 
fits to the broadband data with x^/dof < 1. 

but ion in the variability spectrum. The other would be to 
freeze lnth/l\i— 0 giving a purely thermal spectrum (as in 
AHD13, based on PCA data only). Since this choice does 
not strongly affect the spectral shape below 20 keV both 
variants should give a satisfactory (x^/^ ~1) ht to the data. 
We find that a pure thermal model gives a slightly better 
(x^/^=0'96 compared to 1.05) fit and stick to this model for 
simplicity. We stress that whether the electrons participat¬ 
ing in the fast variability are in fact thermal or hybrid can 
only be determined by higher resolution data (in time and 
energy) at higher energies. Regardless of this choice, how¬ 
ever, the spectrum of the variable component is considerably 
steeper than that of the time-averaged, meaning that even 
if there is non-thermal emission in the spectrum of the vari¬ 


ability, it is not large compared to that of the time-averaged 
spectrum. 

In AHD13 it was found that all the variability spec¬ 
tra of the rising phase of the flaring of XTE J1550-564 had 
little or zero reflection, also this indicating that the vari¬ 
able photons are to a smaller degree intercepted by the disc. 
In our more detailed fit to the co-added frequency resolved 
spectrum we find that reflection, even if low in our best- 
fit (thermal) model, is not compatible with zero. Freezing 
R = 0 only gives an acceptable fit if the absorption line 
is turned into an emission line at 6.45 keV, indicating the 
need for at least a fluorescent iron line even in the variability 
data (regardless of what is assumed regarding the electron 
distribution). Our best-fit thermal model to the spectrum of 
the fast variability has R = 0.24, and requires no absorption 
line. For such low reflection, the ionization parameter can 
not be constrained. We therefore assume the same degree 
of ionization in the reflection of the variability spectrum as 
for the continuum. The best fit model and its components 
is shown in Fig. 4. The parameters are listed in Table 1 
with subscript Var’ as part of the two-component models, 
see next section. (Please note that these values and their 
errors are from the fit to the variability spectrum only. In 
the two-component model fits they are kept frozen.) 

4.3 Two-component comptonization models 

Regardless of the details of the electron distribution and 
exact strength of reflection, the fact that the spectrum of 
the fast variability differs from that of the continuum al¬ 
ready below 20 keV suggests that that there are at least 
two spatially separated electron distributions responsible 
for the comptonized emission of the total spectrum. This 
means that the flow is inhomogeneous and the total con¬ 
tinuum should be the sum of the variable and non vari¬ 
able (on short timescales) part(s). A more physically correct 
broadband spectral model should thus contain at least two 
Comptonization components. This was already briefly dis¬ 
cussed in AHD13 although the data used there (only PCA 
3-20 keV) did not allow for any exact modelling. In DG03 
a two-component model consisting of a hybrid plus an ad- 


© 2013 RAS, MNRAS 000 , ??-?? 



















Pair production and multi-zone comptonization 7 



E [keV] 



Figure 5. Two-component thermal + hybrid Comptonization model to the 0.7-1000 keV data of XTEJ1550-564. One component is 
frozen in shape (but not in strength) to the best fit (fully thermal) model of the spectrum of the 50-100 Hz variability. Left panel: Data 
and model including residuals. ASCA data 0.7-10 keV in black, PCA 3-20 keV in red, HEXTE 20-200 keV in green and OSSE 50-1000 
keV in blue. Right panel: Components of the model. The unscattered blackbody in magenta long dashes, the thermal Comptonization 
component matching the variability spectrum as cyan dot-dashes and its reflection in (shorter) green dot-dashes, Comptonization from 
the hybrid electron distribution in blue short dashes and its Compton reflection as red dots. 




E [keV] 


Figure 6. Two-component thermal + hybrid Comptonization model with full non-thermal injection to the 0.7-1000 keV data of 
XTEJ1550-564. One component is frozen in shape (but not in strength) to the best fit (fully thermal) model of the spectrum of the 
50-100 Hz variability. Left panel: Data and model including residuals. ASCA data 0.7-10 keV in black, PCA 3-20 keV in red, HEXTE 
20-200 keV in green and OSSE 50-1000 keV in blue. Right panel: Components of the model. The unscattered blackbody in magenta 
long dashes, the thermal Comptonization component matching the variability spectrum as cyan dot-dashes and its reflection in (shorter) 
green dot-dashes, Comptonization from a hybrid electron distribution with purely non-thermal injection in blue short dashes and its 
Compton reflection as red dots. The ionization parameters of both reflection components are set equal and determined mostly by the 
stronger hybrid component. 


ditional thermal component was actually found to give the 
best fit to the data (their model HYBTH). The extra com¬ 
ponent in that paper was however ambiguous and as argued 
there may have been only something artificial compensat¬ 
ing for the not-so-correct reflection model. With the help of 
the fast variability we are here able to show that an addi¬ 
tional component should in fact be present in the broadband 
spectrum. 

4 . 3.1 Two component thermal+hyhrid model 

We now fit the spectrum with a model consisting of a 
disc blackbody+ one thermal Comptonization component 


representing the fast variability in addition to the hybrid 
component representing the part of the comptonized spec¬ 
trum that is not variable on fast timescales. We model both 
Comptonized components with eqpair and the total model 
is thus constant*tbabs*gabs(rfxconv*eqpair + KD- 
blur*rfxconv*eqpair). The parameters of the additional 
thermal component are frozen to that of the model for the 
variability spectrum, except for the normalization that is a 
free parameter and allowed larger than or equal to that of 
the variability spectrum which we use as input value. (Keep¬ 
ing it frozen at the same value would be to assume that this 
component is 100 per cent variable). 

To start with, we use the parameter values from the 
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one-component model as input values for the parameters 
of the ’stable’ hybrid component. It turns out that the in¬ 
clusion of the additional component does not make large 
changes to the overall spectrum which is still dominated by 
the hybrid component, the parameters of which are very sim¬ 
ilar in the two models. The normalization of the additional 
thermal component is only slightly higher than the mini¬ 
mum required by the variability spectrum. Besides from it 
being weak, this also indicates that this additional compo¬ 
nent is indeed highly variable. The best-fit two-component 
model is shown in Fig. 5 and its parameters given in Table 
1, column 2. Since the fit is very sensitive to the input val¬ 
ues, we also try to use a higher normalization as input for 
the thermal component. We find that we get still accept¬ 
able {x^/dof < 1.0) but slightly worse fits for the thermal 
component being a factor of ~3 higher (corresponding to it 
being stronger but less variable). The addition of the extra 
component does not alter the limits on Tinj and 7max from 
the one-component model. 


Compton reflected X-rays 



Cool disc 


Hot inner flow, 
hybrid region 


Region of rapid 
variability, thermal 
Comptonization (?) 


Figure 7. Envisaged geometry of the accretion flow with at least 
two Comptonizing regions, one with rapid variability, possibly 
with a purely thermal electron distribution, close to the black 
hole. The rest of the flow which shows strong signs of containing 
both thermal and non-thermal electrons, may in turn be made 
up of several spatially separated regions with different properties 
and electron distributions. 


4-3.2 Two component model with pure nonthermal 
injection 

To try to separate the thermal and non-thermal contribu¬ 
tions, we also model the spectrum with /nth/^h fixed to 1.0 
for the ’stable’ hybrid component. This represents full non- 
thermal injection i.e. that the power supplied in terms of 
injection of high energy electrons is the sole energy sup¬ 
ply to this part of the spectrum. This would be expected 
if this emission is powered by e.g. magnetic flares above 
the disc. (Note that this will not result in a pure non- 
thermal spectrum since the low end will always thermalize. 
Modelling non-thermal injection dominated spectra as pure 
power laws is thus not fully physically motivated even if 
such a model may give a good statistical fit, see e.g. model 
PLTH in GD03 and their discussion). The result is shown 
in Fig. 6 and best-fit parameters for this model are listed 
in Table 1, column 3. In this model the thermal compo¬ 
nent, matching the variability spectrum, and the remaining 
’stable’ non-thermal component are of comparable strength. 
The model requires very high strongly ionized reflection 
and a steep non-thermal injection spectrum, regardless of 
7max, but as before, 7max> 10 overestimates the annihi¬ 
lation line and gives a line flux above the detection limit 
for OSSE. This model does however not give a very good 
fit to the data (x^/dof =1.11, compared to 0.84 for the 
one-component model and 0.82 for the two-component ther- 
mal+hybrid model). It would thus seem that the part of the 
spectrum not variable on short timescales is not compatible 
with being powered solely by non-thermal processes but that 
there are also important contributions from direct heating. 


but are hot when they do causing doppler-smearing of the 
line making it non-detectable. For luminous GBHs the com¬ 
pactness should be high, rather several hundred (see Section 
4) than ~ 10 that is commonly used to avoid the problem 
of creating strong annihilation lines not seen in the data. 
Pair-creation rates are thus very high. With so many soft 
photons present, Gompton cooling is rapid and the pairs 
should lose most of their energy before they annihilate, so 
we do not expect broadening of the line. Since the density 
of pairs is high, created pairs cannot be accreted or escape 
from the source before they annihilate. The only reason¬ 
able explanation is thus a deficit of high energy electrons 
capable of up-scattering photons above the pair production 
threshold. This would require either a very steep injection 
spectrum or a low maximum Lorenz factor of the injected 
electrons, or both. With good quality high energy data, like 
that used here, the the slope of the steady-state electron 
distribution can be rather well constrained and we can de¬ 
termine limits on combinations of Tinj and jmax- The non¬ 
detection of the annihilation line in this data limits jmax 
to < 10 and thus Tinj further to 2.10-2.50 in our best- 
fit models (the one-component and the two-component hy¬ 
brid model). Assuming full non-thermal injection (/nth/G=l, 
two-component non-thermal model) requires steeper injec¬ 
tion Finj= 3.4-4.0, but this model gives a worse fit to the 
data. 

Our results suggest that pair production may not be 
as important in GHBs as previously assumed, at least 
not in the very high state and that preferred acceleration 
mechanisms do not need to produce many electrons above 

'^max — ^^- 


5 DISCUSSION 

5.1 Limits on pair-production 

Despite being predicted by most models, no annihilation 
lines are actually observed in the spectra o f GBHs (except 
for th e claimed detection in IE 1740.7-2942, ISunvaev et al. I 
Il99lh . This can only have three explanations: either no pairs 
are produced, they do not annihilate or they do annihilate 


5.2 Multi-zone comptonization and geometry of 
the flow 

It is by now widely accepted that the origin of the fast vari¬ 
ability in the emission from GBHs in the soft and very high 
state is not the accretion disc but rather the hot flow or 
corona. Our results, here and in AHD13, that the frequency 
resolved spectrum of the 10-50 Hz variability does not con¬ 
tain any disc component confirms this picture. Eurther, the 
fast variability is generally assumed to arise in a region very 
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close to the black hole. Our results, here and in AHD13, 
that the variability spectrum shows less reflection and has 
a harder slope than the continuum indeed suggest that the 
variable region intercepts less soft photons and is located 
far away from the accretion disc. The study of the spectral 
evolution of the variability spectrum versus the continuum 
in AHD13 further showed that the spectrum of the vari¬ 
ability is less sensitive to changes in the accretion rate (and 
thus presumably changes in the inner disc radius) than the 
overall continuum spectrum. Thus, the spectrum of the fast 
variability and its evolution seem to agree with the origin of 
the fast variability being a hot flow close to the black hole 
far away from the disc.The envisaged geometry is shown in 
Fig. [5l For a physical model for how to create the observed 
broad-band variability in the very high state, i nclud i ng th e 
QPO arising in a hot flow, see lingram Don(j(l201ll . l2ni2h . 

Since the variability data covers energies up to 20 keV 
only, the temperature and electron distribution in the vari¬ 
able region can not be constrained but is consistent with 
being purely thermal. To confirm or disprove this we need 
to look at frequency resolved spectra at higher energies. 
Such data are at present not available but may be provided 
by the Hard X-ray Telescope (HXT) aboard the planned 
ASTRO-H mission. The steeper spectrum of the variable 
component, however, already shows that most of the ob¬ 
served non-thermal emission is produced not in the inner¬ 
most variable region, but further out in the flow. It has been 
suggested that the non-thermal emission in accretion flow of 
GBHs are the result of high-energy electrons accelerated in 
magnetic reconnection flares above the disc surface. If this 
is the case then the non-thermal emission arises much closer 
to the disc and should be softer and show stronger signa¬ 
tures of reflection than the variable emission. In our two- 
component models, the ‘stable’ hybrid component is both 
softer and has higher reflection than the variability com¬ 
ponent. We find that this stable component is not likely 
to be powered by purely non-thermal injection but requires 
half of the energy to be supplied as direct heating. This 
is however still compatible with the non-thermal emission 
being produced by magnetic flares above the (truncated) 
disc, since there can also be ’contaminating’ thermal plasma 
from the outer regions of the flow. In EQPAIR thermaliza- 
tion is due to Coloumb interactions only. An even more im¬ 
portant thermalization process, not included in the code, 
is synchrotron self-absorption. The importance of this ef¬ 
fect fo r the emission from compact sources wa s first pointed 
out bv iGhisellini. Guilbert Svensson (llQSSll a n d has been 
further investigated by e .g. IPoutanen VurmI (|2009l j and 
iMalzac BelmontI (|2009l b In the presence of a magnetic 
field the electron distribution may indeed appear thermal 
even if the original acceleration mechanism would have pro¬ 
duced a non-thermal distribution. It is also quite possible 
that the hybrid part of the flow is in fact a combination of 
several regions, each with a different electron distribution. 

We have here investigated the very high state. The sit¬ 
uation is likely to be different in the hard or the more clas¬ 
sical soft state. In the classical soft state (of this and other 
sources) variability is usually strongly suppressed. This is 
consistent with the innermost hot flow collapsing into an 
accretion disc extending all the way to the innermost stable 
orbit. The extra component associated with the fast vari¬ 
ability is thus not expected to be present in the classical 


soft state. In the hard state, variability is stronger and the 
component matching the variability may play a more im¬ 
portant role or even dominate the spectrum. In AHD13 it 
was shown that for the harder of the very high state spectra 
of XTE J1550 the difference between the variability spec¬ 
trum and the total continuum was indeed less than for the 
softe r spectra. Unified spe ctral and timing studies of Gyg X- 
1 bv lYamada et al.l (|2012h have also shown the need for two 
Gomptonized components in the hard state, one of which 
seem to be connected with the fast variability. 


6 CONCLUSIONS 

Galactic black holes are luminous objects and with data cov¬ 
ering most of the soft X-rays from the accretion disc, the 
absolute soft compactness can be calculated, and it is high 
>100. The lack of evidence of any annihilation lines in their 
spectra should thus tell us that perhaps there are less high- 
energy electrons present in the flow than usually assumed. 
We have analysed the broadband 1-1000 keV spectrum of 
the GBH XTE J1550-564. Using a realistic value for the 
compactness we have calculated limits on the high energy 
electron distribution and find that the slope of the OSSE 
data constrains the mean electron energy to be 7 ~ a few, 
requiring either a steep electron spectrum and/or a low max¬ 
imum electron energy. The lack of an observable annihilation 
line favours 7 macc <10 and Finj=2.1-2.5. 

We also model the frequency resolved spectrum of the 
10-50 Hz variability from this source and find that it differs 
from the continuum, not only in that it contains no sign of a 
disc and has very low reflection. It also has a different shape 
of the comptonized part of the spectrum. This means that 
the comptonized flow itself is inhomogeneous and its total 
spectrum must be made up by at least two components, one 
variable on short timescales and one not, where the contin¬ 
uum should be a sum of these. In our interpretation the fast 
variability originates in the innermost parts of the accretion 
flow close to the black hole and far away from the accretion 
disc. This region is spatially and physically separated from 
the origin of the rest of the emission which may in turn be 
a sum of several smaller regions. 

We find that the spectrum of the fast variability is con¬ 
sistent with being fully thermal. This could be conhrmed 
or disproved with access to better quality high-energy data 
from e.g. the HXT onboard ASTRO-H, but our results al¬ 
ready require that the observed non-thermal emission is pro¬ 
duced predominantly further out in the flow. However, the 
time-averaged emission is not consistent with being fully 
non-thermal, indicating that direct heating of the electrons 
is also still important further out in the flow. 
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